Lattice dynamics reveals a local symmetry breaking in the emergent dipole phase of 

PbTe 
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Local symmetry breaking in complex materials is emerging as an important contributor to materi- 
als properties but is inherently difficult to study. Here we follow up an earlier structural observation 
of such a local symmetry broken phase in the technologically important compound PbTe with a 
study of the lattice dynamics using inelastic neutron scattering (INS). We show that the lattice 
dynamics are responsive to the local symmetry broken phase, giving key insights in the behavior 
of PbTe, but also revealing INS as a powerful tool for studying local structure. The new result is 
the observation of the unexpected appearance on warming of a new zone center phonon branch in 
PbTe. In a harmonic solid the number of phonon branches is strictly determined by the contents 
and symmetry of the unit cell. The appearance of the new mode indicates a crossover to a dynamic 
lower symmetry structure with increasing temperature. No structural transition is seen crystal- 
lographically but the appearance of the new mode in inelastic neutron scattering coincides with 
the observation of local Pb off-centering dipoles observed in the local structure. The observation 
resembles relaxor ferroelectricity but since there are no inhomogeneous dopants in pure PbTe this 
anomalous behavior is an intrinsic response of the system. We call such an appearance of dipoles out 
of a non-dipolar ground-state "emphanisis" meaning the appearance out of nothing. It cannot be 
explained within the framework of conventional phase transition theories such as soft-mode theory 
and challenges our basic understanding of the physics of materials. 

PACS numbers: 



I. INTRODUCTION 



The unexpected appearance, on warming, of local Pb 
ofF-centering dipoles was recently reported in PbTe [1]. 
No structural transition is seen in the average rock-salt 
crystal structure but is apparent in the local structure 
on warming above 100 K: the local symmetry is lowered, 
losing its centro-symmetry, on warming. We refer to this 
as "emphanisis" meaning the appearance of something 
from nothing since it is fundamentally different from a 
normal ferroelectric transition where dipoles exist at low 
temperature but become disordered and fluctuating at 
high temperature. Here the dipoles appear on warm- 
ing from a ground-state with no dipoles. In the origi- 
nal study [U the dipoles were deemed to be fluctuating 
but this could not be determined from the experiment it- 
self. A recent inelastic neutron scattering (INS) study [2' 
noted the anharmonicity of certain phonons down to 
low temperature in PbTe. Here we present a detailed 
temperature-dependent INS study of the phonons in the 
temperature range where the local dipoles appear [T| . Be- 
low room temperature our results are in good agreement 
with earlier work on PbTe 3, 4, and we also see the sig- 
nificantly anharmonic signal of the zone center transverse 
optical (TO) mode reported by Delaire et al. P]. How- 



ever, the main result of this work is the characteriza- 
tion of the anharmonic features at ~ 6 meV as a new 
mode with a highly anomalous temperature dependence, 
growing rapidly in spectral weight with increasing tem- 
perature above 100 K, at the expense of the normal TO 
mode that is the incipient ferroelectric mode [4 . The 
new mode and the original TO mode coexist over the 
entire temperature range measured to 600 K, resulting 
in an additional phonon branch in the Brillouin zone in- 
dicative of a broken symmetry, though none is seen in the 
average structure. The new mode is broad with a short 
lifetime, but dispersive. It also hardens with increasing 
temperature. Since it appears over the same temperature 
range where local dipoles appear in the structure [T] we 
associate its appearance with the appearance of these ob- 
jects that break the local, though not the average, centro- 
symmctric symmetry. Such behavior, observed in a pure 
binary alloy system, has not been described before and 
challenges our current understanding of the physics of 
materials. 



II. RESULTS 

We first consider the lattice dynamics at room tem- 
perature and below, and compare our measurements to 
earlier results in the literature [5H1] to establish the qual- 
ity of our sample and data. Fig. [I]ja)-(c) shows repre- 
sentative room temperature INS constant Q scans, at 
three points in the Brillouin zone, collected on the HB3 
triple axis spectrometer at Oak Ridge National Labora- 
tory (ORNL). Details of the data collection are described 
in the Materials and Methods section below. As illus- 
trated in the inset to Fig.[T] the points in reciprocal space 
are (a) (033), which is a zone boundary iiT-point ^ in a 
position such that the instrument is sensitive to trans- 
verse polarized modes [6] (b) (2.5 2.5 2.5), which is a 
zone boundary L-point sensitive to longitudinal modes, 
and (c) (133), which is a Brillouin zone center (F-point) 
where both longitudinal and transverse modes can be 
measured. Certain phonon branches are expected based 
on the earlier work of Cochran et. al [3], the energy 
transfers of which are indicated by red dashed lines in 
Fig. [11 There is good semi-quantitative agreement with 
the earlier results. Fig.[I](d)-(f) shows a false-color plot of 
the scattered intensity as a function of temperature and 
energy transfer, AE = fku, for the same three Q-points. 
The ii'-point modes in Fig. flld) do not change energy 
at all vs. temperature, neither do they broaden signifi- 
cantly. This is expected in a harmonic system, but even 
in a quasi-harmonic approximation where lattice expan- 
sion is taken into account it is common to see softening 
with increasing temperature [7 . The LA/LO doublet at 
the L-point shown in Fig. W[c) does appear to broaden 
with increasing temperature though the peak contains 
two components that are both softening with increasing 
temperature. Fits to the line-shapes suggest that there 
is little broadening of the individual lines. Qualitatively, 
the softening of these modes may be explained by a quasi- 
harmonic analysis. 

We now consider the zone-center TO mode in Fig. IlTc) 
and (f). This mode is the soft-mode that indicates an in- 
cipient ferroelectric phase transition in PbTe [1]. The 
mode softens with decreasing temperature (the opposite 
of the L-point modes and opposite to early density func- 
tional theory (DFT) calculations [8]) but never reaches 
zero frequency. A softening to zero frequency would in- 
dicate the soft-mode transition temperature ^ and a 
structural phase transition. It is immediately apparent 
that the mode is anomalously broad as pointed out re- 
cently [T. The zone edge modes in Fig.fTJa) and (b) have 
a FWHM of around 1.5 meV which is close to the calcu- 
lated energy resolution of the instrument, as described in 
Materials and Methods below. However, the high-energy 
tail on the zone center TO mode is of order ~ 5 meV. 
The temperature dependence of this scattering feature is 
shown in a false-color plot in Fig. flVf ) . This feature in 
the scattering is seen to broaden dramatically with in- 
creasing temperature and is highly asymmetric and non 
Gaussian [2]. 



To study the temperature dependence of this mode in 
greater detail we have fit curves to the spectra at each 
temperature. We focus initially on the low-temperature 
region, T < 300 K. We would like to see if our data 
reproduce the earlier results of Alperin et al. |4j. Rea- 
sonable fits could be made to the scattering features by 
fitting a Lorentzian peak on a linear sloping background, 
and from this we extracted the mode frequency. Three 
of these fits are shown in FigsJl|g)-(i). The temperature 
dependence is shown in Fig. fTm. For a mean-field, sec- 
ond order phase transition, the soft mode theory predicts 
that the frequency of the soft mode approaches zero as 
(T - Tcf-^ 0, so a plot of w2 vs. T wih be linear and 
intercept the abscissa at transition temperature, Tc. It 
is clear that the zone center mode softens on cooling ac- 
cording to mean-field soft mode theory but still has a 
positive frequency at T = 0: PbTe is an incipient ferro- 
electric but, in the absence of alloying with Ge [TO^, does 
not undergo a ferroelectric phase transition on cooling. 
There is excellent agreement with the earlier work [3]. 
Extrapolating the line to the abscissa we obtain a pu- 
tative phase transition temperature of -151 K. This is a 
little more negative in temperature than Alperin et al. @] 
(-135 K) but in good general agreement and establishes 
the purity of our crystal. 

Having established that data from our sample are in 
good agreement with earlier work we move to the main 
result of the current study. In Fig. [2] we show the temper- 
ature dependence of the zone center TO mode feature at 
and above room temperature. Over a narrow tempera- 
ture range beginning above ~ 100 K a new peak grows up 
rapidly but smoothly from the high energy tail of the TO 
mode, gaining spectral weight at the expense of the TO 
mode, but coexisting with it. Above room temperature 
the new peak becomes very strong, dominating the orig- 
inal TO feature. It indicates the presence of a new mode 
in the system that appears on increasing temperature, 
which is highly anomalous. 

The neutron scattering dynamic structure factor, 
S{Q,ll)), is the normalized inelastic scattering intensity 
and is related to the imaginary-part of the dynamic sus- 
ceptibility according to [6 
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For a particular phonon mode in the sth branch at posi- 
tion Q = G—q'vn reciprocal space, where G is a reciprocal 
lattice vector and g is a vector that lies within a single 
Brillouin zone, x"((3,w) is proportional to, \F{Q)\'^ /uj^g, 
where F{Q) is the phonon dynamic structure factor. It is 
proportional to the Bragg scattering structure factor for 
zone center acoustic modes, but in general depends on 
the polarization of the vibrational mode and the atomic 
displacement vectors of the eigenfunction. 

To obtain a quantity proportional to \" {Q, w), we mul- 
tiply the measured intensity by the Bose-Einstein factor 

C — (1~ e ''B^ ) . This removes the temperature depen- 
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FIG. 1: (Color) INS spectra of a PbTe single crystal at 288 K at different Q points in reciprocal space: (a) (0 3 3) zone 
edge, (b) (2.5 2.5 2.5) zone edge, and (c) (1 3 3) zone center. The red dashed lines indicate the phonon energies measured 
by Cochran et. al [3] for comparison. The inset in (c) is a schematic of the plane of reciprocal space that the triple axis 
spectrometer was set to work in. The grey dots are the reciprocal lattice points and the lines are Brillouin zone boundaries. 
The red, pink and blue dots are the Q points for the data shown in (a) (A'-point), (b) {L point) and (c) (F-point), respectively, 
(d)- (f) show false color plots of the intensity (color axis) vs. energy transfer, AE — huP' , and temperature of the same modes 
shown in the panels above. (g)- (i) show representative fits of Lorentzian peaks (in red) to the TO mode at the [133] zone 
center spectra: (g) 2 K, (h) 134 K and (i) 231 K. The vertical dashed lines are the extracted mode energies, cj. (j) uP vs. T 
for all the temperatures fit in the temperature range below 300 K with a straight line fit to the data shown in red. 
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FIG. 2: (Color) Plots of the INS spectrum at Q = (133). 
From the bottom (pale blue) the temperatures shown are 
28 K, 99 K, 203 K, 299 K, 410 K, 486 K, and 601 K (topmost 
curve). Datasets are offset with respect to each other by 75 
intensity units for clarity. The two arrows are guides to the 
eye, indicating the positions of the TO and the new mode at 
the highest temperature. Both modes soften on cooling. 



dent effects of phonon occupancy and allows us to track 
the changes in the phonon energy spectrum directly . 
The thus normalized datasets from Q = [331] were then 
fit with 2 damped harmonic oscillator functions to ex- 
tract the different phonon dynamic susceptibilities. An 
additional Gaussian function was introduced to fit the 
intensity from the tail of the Bragg peak. Details of the 
fitting are in the Materials and Methods section. Rep- 
resentative examples of the fits to the data at various 
temperatures are shown in Fig. p[a)-(f). The growth of 
the susceptibility of the new mode, at the expense of the 
TO mode, with increasing temperature is clearly appar- 
ent, as evident in Fig. |3[g). As well as growing, the new 
mode sharpens and hardens with increasing temperature. 
Rapid growth of the new mode spectral weight begins at 
around 100 K. The behavior is reminiscent of the power 
law growth of an order parameter in a Ginsberg-Landau 
theory of a conventional second order phase transition [5^ 
except in that case the order parameter grows upon cool- 
ing through the phase transition, not on warming as here. 
The temperature dependence of the mode susceptibil- 
ity is reminiscent of the growth in the amplitude of the 
off-center Pb displacements observed in the atomic pair 
distribution function (PDF) [I], as shown in the inset. 
We have investigated the dispersion of the new mode 
in the [/ill] direction away from the zone center using 
the ARCS chopper spectrometer at ORNL. As is evident 



from Fig. [4J the mode clearly disperses to higher energy 
transfer moving away from the zone center, similar to the 
TO mode, and hardens with increasing temperature. 

The observation of additional modes implies a symme- 
try breaking in the system. The most trivial possibility is 
the appearance of a localized mode due to the presence of 
point defects in the material. Such modes are extremely 
weak and non-dispersive but can become apparent, even 
for modest defect densities, when they interact with dis- 
persing crystal modes [11] . This explanation can be ruled 
out by the temperature dependence. 100 K is too low a 
temperature for any significant thermally activated de- 
fect formation and it is difficult to explain the power- 
law behavior of the susceptibility in this scenario. Also, 
the new mode is dispersive arguing against this explana- 
tion. Similar arguments argue against another possibility 
which is the formation of thermally induced intrinsic lo- 
cal modes [121. These are localized phonon modes that 
form due to non-linear effects in a strongly driven system 
of oscillators [13] and have been postulated to form due to 
thermal excitation in soft ionic systems [12]. An anhar- 
monic coupling between the LA and TO modes has also 
been postulated to explain the breadth of the TO mode 
in PbTe [5], but this does not explain the characteris- 
tic temperature dependence, and the clear observation of 
two distinct coexisting components in the vicinity of the 
TO mode. Also the LA mode is at zero frequency and far 
from the TO mode at the zone center where our data were 
measured. The most striking correspondence of the new 
mode is with the temperature dependent appearance of 
Pb off-centering structural dipoles in the local structure 
reported from PDF measurements [J. This observation 
was rationalized on thermodynamic grounds as the en- 
tropically stabilized appearance of a paraelectric phase 
above an undistorted non-ferroelectric ground-state. In 
this picture, the emergent dipoles break the local symme- 
try by removing the center of symmetry on the Pb site, 
even though there is no change in the long-range sym- 
metry. This will result in short and long Pb-Te bonds, 
changing mode frequencies associated with the TO vibra- 
tions which are the intra-unit cell anti-phase vibrations 
of Pb and Te, just as observed here. If each off-center 
Pb atom is fluctuating independently, this would result 
in a flat non-dispersing mode. The observation of dis- 
persion implies that off-centered Pb ion displacements 
are correlated over some range of space, resulting in po- 
lar nanodomains similar to those observed in relaxor fer- 
roelectrics below the Burns temperature [T4|. The new 
mode seen here is highly reminiscent of modes appear- 
ing in INS measurements of relaxor ferroelectrics that 
were attributed to polar nanodomains |15| . However, 
the crucial distinction here is that PbTe is a pure ma- 
terial. In relaxors the polar nanodomains are stabilized 
by nanoscale chemical disorder whereas here they must 
be intrinsic. A physical explanation may lie in the par- 
tially localized optical phonon modes recently postulated 
byDFT and molecular modeling [16j that result from par- 
tial covalency of the Pb and Te that competes with the 
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FIG. 3: (Color) (a)- (f) : Dynamic susceptibility obtained from the measured INS spectra (symbols) at representative 
temperatures with fits to the data shown as a pink line. There are 2 damped harmonic oscillator components, one for the TO 
mode (red) and one for the new mode (green). A Gaussian function was added to fit the intensity from the tail of the Bragg 
peak (orange). The horizontal blue line is a constant background component, (a) 2 K, (b) 98 K, (c) 203 K, (d) 299 K, (e) 
482 K and (f) 601 K. (g) shows the temperature dependence of the integrated weight in the new mode from the fits (black 
symbols). The inset shows the Pb off-centering displacement measured from the PDF measurement on PbTe reported in jlj 
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FIG. 4: (Color) Neutron scattering intensity vs. energy transfer from single crystal PbTe measured on the ARCS spectrometer. 
Curves in each panel are cuts along [/ill] for 5 (zone center)< h < 5.2 showing the dispersion of the modes, (a) 100 K, (b) 
200 K, (c) 300 K, (d) 400 K, (e) 500 K and (f) 600 K. At high temperature the intensity is almost exclusively coming from the 
new mode which is seen clearly dispersing to higher energy with increasing h. 



electrostatic Madelung potential which prefers the high 
symmetry rock-salt structure. More work on the tem- 
perature dependence and the dispersion of these modes 
is needed to establish the true origin of this behavior. 

We would like to establish if there is a static compo- 
nent to the short-range ordered nanoscale distortions by 
searching for diffuse scattering around the zone center 
in the elastic channel (w = 0). Indeed, we find there is 
a small increase in scattering at the base of the (111) 
Bragg peak at high temperature, which could be elas- 
tic diffuse scattering, as evident by comparing Figs, pfa) 
and (b) . The Bragg peak is evident as the strong scatter- 
ing signal in dark red centered at (1.0,1.0,1.0). Another 
strong feature at (1.25,1.3,1.3) is a spurious feature of 
unknown origin. Powder diffraction lines are also evi- 
dent coming from the aluminum cryostat windows, and 
for the line going through the Bragg peak, possibly some 
crystal mosaicity. The spurious features disappear (or 
rather appear with negative differential intensities) in 
the difference plot in Fig. Isjc) which shows the differ- 
ence between the intensities at high, (b), and low, (a), 
temerature. We see some diffuse scattering intensity ap- 
pearing at high temperature in this elastic measurement, 
suggesting that the high temperature broken symmetry 
phase has a static component. This extra diffuse scatter- 
ing component shows up as positive intensity around the 
(111) point in Fig. [5|c) whereas the Bragg features have 
a negative difference (dark blue) due to Debye- Waller ef- 
fects. However, the most likely origin of the observed 
increase in elastic diffuse intensity is the finite energy 
resolution of our measurement resulting in a temperature 
dependent signal coming from the increase in population 
of the acoustic phonons that lie within the energy reso- 
lution window. Fig. p[d) and (e) show scans that have 
been approximately corrected for phonon effects by di- 
viding the region around the Bragg peak by the Bose 
factor appropriate for the phonon modes in the resolu- 
tion window. This correction accounts for all of the ob- 
served diffuse scattering at high temperature in Fig. [5[c) 
(see Fig. IsFf)) suggesting that there is no elastic diffuse 
scattering signal from the dipoles and therefore they have 
no static component but are completely dynamic in na- 
ture, again consistent with the DFT prediction [TB] 

The chemically simple binary PbTe compound contin- 
ues to turn up surprises that challenge our understand- 
ing of condensed matter systems. We know of no other 
example of the appearance in a pure material of a new 
phonon mode with increasing temperature, and the char- 
acteristic temperature dependence of its dynamic suscep- 
tibility, argues strongly that this is an intrinsic response 
of the system. The inelastic scattering is sensitive to 
the local symmetry which is broken at high tempera- 
ture where the new mode is seen to coexist over a wide 
range of temperature with the original TO mode of the 
undistorted structure. A picture emerges of nanoscale 
regions where the local dipoles appear at high temper- 
ature and are correlated and fluctuating. This lowering 
of local symmetry with rising temperature, betrayed by 



the emerging modes seen here and in the PDF [T], may 
explain the long known anomalous temperature depen- 
dence of the semiconducting energy gap in PbTe [121 [H] . 
Unlike in conventional semiconductors where the energy 
gap is known to increase with falling temperature due 
to the reduction of the intensity of thermal vibrations, 
that of PbTe increases with increasing temperature up 
to as high as 600 K. It is this anomalous dependence 
that permits PbTe to exhibit delayed cross gap carrier 
excitations thereby sustaining a very high thermoelectric 
power factor at high temperatures. Coupled with the in- 
creased anharmonicity and phonon scattering from the 
same origin resulting in a very low thermal conductivity, 
this may help to explain the one two punch that propels 
PbTe to the top of its thermoelectric class TE] . 



III. MATERIALS AND METHODS 

A. Sample preparation 

Stoichiometric amounts of Pb (Rotometals, at 99.9% 
purity) and Te (Plasmaterials, at 99.999% purity) were 
flame sealed in an evacuated (< 10~^ mTorr) fused silica 
tube. The mixture was heated to 1050 C at a rate of 
100 C/h for 4 h and cooled down to room temperature 
at a rate of 20 C/h. For the single crystal growth, ~ 20 g 
of PbTe were remelted using the Bridgman method. A 
13 mm fused silica tube was loaded with PbTe and low- 
ered at a rate of ^ 2.6 mm/h through a single zone verti- 
cal furnace that was set at 1050°C. A cylindrical shaped 
crystal of 13 mm in diameter and ~ 60 mm in height was 
used for the inelastic neutron scattering experiment. The 
quality of the single crystal was evaluated using neutron 
Laue diffraction at SNS. 



B. Triple Axis Inelastic Neutron Scattering 
Experiments 

The inelastic neutron experiments were performed at 
the IIB3 triple axis instrument at the High Flux Iso- 
topie Reactor (HFIR) reactor located at Oak Ridge Na- 
tional Laboratory. All the experiments were performed 
with a fixed final energy, Ef— 14.7 nieV using Pyrolytic 
Graphite (PC) as a monochromator and analyzer and a 
PG filter after the sample to eliminate higher order con- 
tamination of the scattered beam. The neutron beam 
coUimation was 48-20-40-90 before the monochromator, 
sample, analyzer, and detector, respectively. This yielded 
an energy resolution of 1.06 meV full width half max- 
imum (FWHM) for zero energy transfer, increasing to 
1.62 meV FWHM at w = 12.0 meV energy transfer. The 
sample was mounted to the cold-tip of a closed cycle He- 
4 based refrigerator (CCR). Data scans were made as a 
function of energy transfer at constant Q points in recip- 
rocal space at 35 different temperatures between 2.5 K 
and 600 K. 



A 

i.5|_rt''- 

1.4 



1.3 

^1 2 

1.1 

1.0 











9-r'v r -J 



0.8 



1.5 







1.1 u_ 

1.0 

0.9 tii !! 

0.8 










0.8 1.0 1.2 1.4 O.i 
hOO 



1.0 1.2 1.4 



JIO 




■"' ■yi.'l".-'"' 







^^^^'^ ■; 



I 



20 
10 § 

< 





_>> 
I- -I-10S 

20 



■ ti I " B m 



0.8 1.0 1.2 1.4 

hOO 



I 

I 

I 



20 

10 § 
xi 

^ 

>. 

'w 
-10| 

c 

-20 



FIG. 5: (Color) Elastic diffuse scattering around the [hOO],[Okl] plane of reciprocal space measured using the Triple axis 
spectrometer: (a) 16 K, (b) 482 K (c) difference between the scans in (a) and (b): 1(482) — /(16). (d) and (e), as (a) 
and (b) but the intensities have been corrected to account for the change in phonon occupation of the acoustic modes that are 
within the energy resolution of the measurement, (f) The difference between the scans in (d) and (e) 



C. Triple Axis Inelastic Neutron Scattering Data 
Curve Fitting 

Here we describe the method used to extract the dy- 
namic structure factors of the modes at the (133) zone 
center from the TAS data. First the data are normal- 
ized to obtain the dynamic susceptibihty, ^"(Qj'^)? ^ 
described in the main paper. Based on earher work we 
expect two zone-center modes to be present at finite en- 
ergy transfer (the acoustic modes are hidden under the 
elastic line at w = 0): the TO mode at around 4 meV 
and the LO mode around 14 meV. As the LO mode does 
not overlap with either the TO or the new mode, we do 
not take this into account in the fit, which only includes 
the region from 1-12 meV. The TO mode is then fitted 
with a function representing a damped harmonic oscilla- 
tor [H], with an additional constant background in the fit. 
Furthermore, a Gaussian function is added to take care 
of the tail of the Bragg peak which is present even after 
the normalization with the Bose factor. We then add 
an additional damped harmonic oscillator to account for 
the new mode. The spectra were then fit as a function 
of temperature sequentially with the starting model for 
each fit using the parameters obtained from the previous 
temperature point. Certain constraints were introduced 
in order to aid convergence, namely, the widths of TO 
and new mode curve and the background level were al- 
lowed to change by only ±20%, between subsequent runs. 



D. Triple Axis Elastic Diffuse Scattering 
Measurements 



Elastic scans were also taken on a grid of points in the 
Brillouin zone At 16 K and 482 K to search for diffuse 
scattering from static disorder. The range 0.8 < h < 1.5 
was scanned in the [hOO] direction and 0.8 < k < 1.5 
in the [Okk] direction to capture the region around the 
(111) reciprocal lattice point. This is a zone center and 
the PbTe Bragg peak is evident in Fig. Isla) at the (111) 
point. There is an additional Bragg feature at around 
(1.25 1.3 1.3). We are not sure of the origin of this fea- 
ture and think it is a spurious peak. It does not show any 
interesting temperature dependence and we neglect it in 
the analysis. In addition there are a number of powder 
Bragg lines evident. One goes through the (111) point 
and suggests some sample mosaicity. The other powder 
lines probably come from the aluminum sample environ- 
ment. 

In Fig. p[d) and (e) we have applied a simple correc- 
tion to account for the increased occupancy of acous- 
tic phonons that lie within the energy resolution of the 
measurement. Based on the dispersion curves measured 
by Cochran et. al, we estimate the range of Q that 
is affected by phonons within the energy resolution (of 
around 1 meV). We assume that the average phonon 



energy in that region is ~ 0.5 meV and so correct the 
intensities in this range of Q by the Bose factor appro- 
priate for a mode of that energy and the temperature 
in question. This is a highly simplified correction and is 
intended to estimate the scale of intensity coming from 
the increase mode population at high temperature rather 
than being a highly accurate correction. We see that this 
simple correction accounts for essentially 100% of the dif- 
fuse scattering seen at higher temperature in the elastic 
channel. 



E. Chopper Spectrometer Inelastic Neutron 
Scattering Experiments 

Inelastic neutron scattering measurements were also 
performed using the ARCS direct geometry time-of-fiight 
chopper spectrometer at the Spallation Neutron Source 
(SNS) at the Oak Ridge National Laboratory. Single 
crystal measurements were performed with an incident 
energy of Ei = 30 meV with the sample mounted approx- 
imately in the (hkO) scattering plane of the instrument. 
The sample was mounted to the cold-tip of a closed cy- 
cle He-4 based refrigerator (CCR). Data were acquired by 
rotating the sample about the vertical axis by 40 degrees, 
and collecting data in 1.875 degree increments. Empty 
sample can measurements were performed for all data ac- 
quired at ARCS and subtracted from the data shown in 
the manuscript. 

The resulting large 5-dimensional data-sets show the 
scattering intensity throughout much of {Q,ijj) space. 
The data are projected into a 4D space and 3D slices 
and 2D cuts are made to extract features of interest using 
the DCS-mslice program within the DAVE [T^ software 
package. 
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